Magnetoresistance of magnetic multilayers in the CPP mode: 
evidence for non-local scattering 
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We have carried out measurements of the magnetoresistance MR(H) in the CPP (Current Per- 
pendicular to the Plane) mode for two types of magnetic multilayers which have different layer 
ordering. The series resistor model predicts that CPP MR(H) is independent of the ordering of the 
layers. Nevertheless, the measured MR(H) curves were found to be completely different for the fol- 
lowing two configurations: [Co(10A)/Cu(200A)/Co(60A)/Cu(200A)]n and [Co(10A)/Cu(200A)]n 
[Co(60A)/Cu(200A)]n showing that the above model is incorrect. We have carried out a calculation 
showing that these results can be explained quantitatively in terms of the non-local character of the 
electron scattering, without the need to invoke spin-flip scattering or a short spin diffusion length. 



PACS numbers: 75.70.Pa, 75.70.-i, 73.40.-c 

Since the discovery a decade ago of the giant mag- 
netoresistance exhibited by magnetic multilayers, inter- 
est in this phenomenon has not abated ||l|. Recent re- 
search has focused on the magnetoresistance MR(H) in 
the CPP mode (current perpendicular to the plane of 
the layers) ^-^. Measurements of MR(H) are techni- 
cally more difficult in the CPP mode than in the CIP 
mode (current in plane). However, there are advantages 
to the MR(H) data in the CPP mode. For example, it 
has been shown [Q||] that experimental values of MR(H) 
in the CPP mode can shed light on the spin diffusion 
length. Here we present evidence for the importance of 
MR(H) measurements in the CPP mode for determining 
the role of non-local electron scattering in the giant mag- 
netoresistance (GMR). We shall show that because of the 
long electron mean free path, non-local scattering makes 
the series resistor model inappropriate. 

As is well known, the GMR occurs in magnetic mul- 
tilayers because the spin-up electrons and the spin-down 
electrons have different scattering rates. If the electron 
does not flip its spin upon scattering, then the spin-up 
and spin-down electrons constitute two separate currents, 
with different resistivities, as if flowing in two parallel 
wires. In the CPP mode, the resistances of the different 
layers add in series |Q,0|g]. Therefore, it would seem that 
two magnetic multilayers that differ only in the ordering 
of the layers would yield identical results for MR(H) in 
the CPP mode. 

To test this idea, Pratt and co-workers at Michi- 
gan State University (Chiang et al. |g[) measured CPP 
MR(H) for the two configurations [Py/Cu/Co/Cu]N and 
[Py/Cu]N[Co/Cu]N (denoted as 'interleaved' and 'sepa- 
rated' configurations, respectively), where Py is Nig4Fei6. 
Although the expectation was that identical MR(H) 
curves would be obtained for the interleaved and the 



separated configurations, these workers found that the 
resulting two MR(H) curves were completely different. 
Chiang et al. [g[ attributed their results to the short spin 
diffusion length in Py. They had previously analyzed re- 
sistivity data within the framework of Valet-Fert theory 
|0,|| and obtained |10| for Py a spin diffusion length of 
only 55 A, thus implying significant mixing between the 
spin-up and spin-down electron currents. Chiang et al. 
proposed that this spin-flipping was responsible for the 
different CPP MR(H) curves they observed for the sepa- 
rated and interleaved configurations. 

We have investigated these ideas by measuring MR(H) 
for multilayers whose magnetic layers do not exhibit 
a short spin diffusion length. For the different mag- 
netic layers, we used Co of two different thicknesses, 
since Co is known |f^Jf3] to have a long spin diffu- 
sion length. Measurements were carried out of CPP 
MR(H) for [Co(I0A)/Cu(200A)/Co(60A)/Cu(200A)]n 
and [Co(I0A)/Cu(200A)]n[Co(60A)/Cu(200A)]n for 
N = 4, 6, 8. The thickness (200 A) of the non-magnetic 
layers was chosen to be large enough to ensure complete 
magnetic decoupling between the ferromagnetic layers. 
In spite of the fact that the interleaved and separated 
configurations differ only in the ordering of the layers, the 
measured MR(H) curves were found to be very different 
for the two different configurations. We shall show that 
these results can be explained quantitatively in terms of 
non-local electron scattering. 

The multilayers were grown in our VG-80M MBE facil- 
ity which has base pressure of typically 4x10^^^ mbar. 
Our CPP measurements used the superconducting Nb 
electrode technique, as developed by Pratt et al. |g] . The 
superconducting equipotential [0,0 ensures that the cur- 
rent is perpendicular to the layers. We used a SQUID- 
based current comparator, working at 0.1% precision to 
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FIG. 1. Magnetic field dependence of tfie magnetoresis- 
tance MR(H) for tlie interleaved (squares) and separated (cir- 
cles) multilayers containing Co (10 A) and Co (60A) as the two 
magnetic metals, for the indicated number of repeats. 



measure changes in the sample resistance of order 10 pJ7. 
To avoid driving the Nb normal, the CPP measurements 
were performed at 4.2 K in magnetic fields below 3 kOe. 
Consistency between the interleaved and separated sam- 
ples was enhanced by growing the two configurations dur- 
ing the same run for each value of N. 

The magnetoresistance was measured in the CPP 
mode for the two configurations: interleaved and sepa- 
rated. The measured curves for MR(H) are presented for 
three values of N in Figs, la-lc. The squares represent 
the MR(H) data in the interleaved configuration whereas 
the circles give the data in the separated configuration. 
For each sample, the saturation magnetic field was about 
2 kOe. 

There are several characteristic features of these data, 
all of which can be explained in terms of non-local elec- 
tron scattering, (i) The most important feature is surely 
the striking difference between the MR(H) curves for the 
two configurations, both in shape and in magnitude, 
(ii) For each N, the maximum value of MR(H) is larger 
for the interleaved configuration, (iii) The MR(H) curve 
for the interleaved configuration exhibits a single peak, 
whereas for the separated configuration, MR(H) is the 



superposition of two peaks, with the second being much 
broader and less delineated than the first. Another in- 
teresting feature of the data, not displayed in Fig. 1, is 
that the saturated resistance itself is always greater for 
the interleaved configuration. 

To ensure that the differing results for MR(H) for the 
two configurations are not due to differences in their mag- 
netic properties, the magnetization as a function of field 
was measured for each sample. We found that the two 
configurations yield the same magnetization. This con- 
firms that the magnetic layers are uncoupled and become 
magnetized independently. At low fields, the magneti- 
zation curves are dominated by the contribution of the 
thicker Co layers. After the thicker Co layers reach sat- 
uration, the magnetization continues to increase as the 
thinner Co layers approach saturation. The magnitudes 
of the saturation fields for the two thicknesses of Co lay- 
ers correspond closely to the saturation fields of MR(H) . 

Kinetic theory arguments show that the electron mean 
free path is far longer than the thicknesses of the mag- 
netic layers (10 A and 60 A). Therefore, the potential 
"felt" by the electron is the combined potential of a 
neighboring pair of magnetic layers. This may be termed 
" non-local" electron scattering in the sense that one can- 
not speak of the resistivity of a single Co layer. Rather, 
the resistivity is determined by a property of pairs of 
neighboring layers. Gittleman et al. [13[ have shown that 
for such a case, the contribution of the spin-direction- 
dependent resistivity depends on the cosine of the angle 
9ij between the moments of neighboring magnetic layers, 
i and j. This is the key to understanding the data. 

Because the mean free path is larger than the layer 
thicknesses, it is necessary to carry out a full band struc- 
ture calculation to calculate properly the resistivity and 
magnetoresistance. However, one can understand the ba- 
sic physics with a simple phenomenological model. 

For the interleaved configuration, the neighboring 
magnetic layers are different, and hence the maximum 
angle 9ij is large, whereas for the separated configura- 
tion, the neighboring magnetic layers are the same (ex- 
cept for one boundary layer), and hence the maximum 
angle 6*^ is small. Therefore, there is no reason to expect 
MR(H) to be the same for the two configurations. This 
explains the first feature of the data mentioned above. 

From the above considerations, it also immediately fol- 
lows that MR(H) will be larger for interleaved multilayers 
than for separated multilayers, because the angle 9ij is 
larger for the former configuration. This explains the sec- 
ond feature of the data mentioned above. This has been 
confirmed by measurements of the GMR as a function of 
the number of bilayers. A Fuchs-Sondheimer analysis of 
these data shows that the mean free path in sputtered 
^ and MBE ^ samples is about 500 A and 700 A, 
respectively. 

For the interleaved configuration, there is only one an- 
gle 9ij that is relevant, namely, the angle between the mo- 



merits of the different (lOA and 60A) neighboring mag- 
netic layers. Therefore, there will be only one peak, as 
the angle Oij becomes progressively larger, passes through 
a maximum at the saturation field of the Co (60A) layer 
and then becomes smaller as the Co (lOA) layer also 
saturates. By contrast, for the separated configuration, 
there are two angles Oij that are relevant, namely, the 
angle between neighboring moments for each set of lay- 
ers (the 10 A set and the 60 A set). As each angle Oij 
passes through its maximum, a peak will be obtained 
for MR(H), leading to two overlapping peaks, with each 
maximum occurring at a different value of the magnetic 
field, corresponding roughly to the coercive field of each 
type of magnetic layer. This explains the third feature 
of the data mentioned above. 

These ideas can be made quantitative. If the spin dif- 
fusion length is very long, it is known [y that a sim- 
ple expression is obtained for MR(H). According to the 
phenomenological theory of Wiser ||l^ , for the geometry 
under consideration here and assuming a very long spin 
diffusion length, the magnetoresistance due to an ij-pair 
of neighboring magnetic layers is: 



MR,j{H) = cy (1 - cos 0,,{H)Y 



(1) 



The spin diffusion length of Co has been measured 
yielding values of 450 A ^j^ and 1000 A §. These 
values is very much larger than the thickness of the Co 
layers, and so one may safely employ the expression for 
MR(H) given in (1). 

For our samples, there are three parameters Cij corre- 
sponding to the three different types of neighboring pairs 
of magnetic layers: i = j = 1; i = j = 2; i = 1, j = 2, 
where 1 refers to Co (60 A) layers and 2 refers to Co (lOA) 
layers. The interleaved configuration contains only type 
i = 1, j = 2 neighbors, whereas the separated configura- 
tion contains all three types. For a sample containing N 
repeats, the separated configuration consists of N-1 pairs 
of type i = j = 1 neighbors, followed by one pair of type 
i = 1, j = 2 neighbors (the boundary layer), followed by 
N-1 pairs of type i = j = 2 neighbors. 

First consider the interleaved configuration. The sat- 
uration magnetic field Hgi of the thicker Co layers is 
smaller than Hs2 of the thinner Co layers. Thus, as the 
magnetic field is increased, the angle Oi^2 increases, since 
the thicker Co layers are reversing their direction of mag- 
netization faster than the thinner Co layers. According 
to Eq. (1), increasing the angle 9i^2 implies an increase in 
MR(H). When the magnetic field reaches Hgi, the angle 
9i2 reaches its maximum value, and begins to decrease as 
the thinner Co layers continue to reverse their direction 
of magnetization while the thicker Co layers have already 
reached saturation. According to Eq. (1), decreasing 0i^2 
leads to a decrease in MR. Finally, when the field reaches 
Hs2, the angle 61^2 is again zero, and MR vanishes. Thus, 
we expect - and find - a single peak for MR(H) for the 
interleaved configuration. 
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FIG. 2. Comparison between the calculated curves (solid 
lines) and the data points (squares) for MR(H) for the inter- 
leaved configuration for the indicated number of repeats. 



The field dependence of Oij is determined as follows. 
The magnetization increases linearly with field (except 
near saturation, where it increases more slowly). Since 
the magnetization is proportional to the cosine of the an- 
gle between the magnetic moment and the field, it follows 
that cos^i and cosOj are each linear in the field. Equation 
(1) contains cosOij — cos{Oi - Oj). Expanding the cosine 
readily gives the required field dependence. 

The calculated results |0 for the interleaved configu- 
ration are given by the curves in Fig. 2. For each value 
of N, the parameter ci.2 was determined by fitting to 
the MR(H) data. The agreement between the calculated 
curves and the data is evident from the figure. 

We now consider the separated configuration. If the 
Co layers were ideal single-domain structures, then the 
magnetic moment of each Co layer would react identically 
to the magnetic field and the angles Oi^i and 02,2 would 
both be zero at all fields. However, because of the pres- 
ence of domains and of structural imperfections in the Co 
layers, each layer reverses its magnetization at a some- 
what different rate. As a result, the angles Oi^i and 6*2.2 
become non-zero as the field is increased, pass through a 
maximum at the coercive field, and then decrease to zero 
as saturation is approached. 

We assumed a simple parabolic form for each of the two 
angles. The maximum value of each parabola, Omax,i,i 
and Omax,2,2, cannot be determined by fitting to the data 
for the following reason. Because these angles are small, 
Eq. (1) can be expanded to yield 
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FIG. 3. Comparison between the calculated curves (solid 
lines) and the data points (circles) for MR(H) for the sepa- 
rated configuration for the indicated number of repeats. 
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and this combination of cu and Omax,ii serves as a single 
fitting parameter. Nevertheless, some numerical tests we 
have carried out suggest that both 0max,i,i and 0max,2,2 
lie in the range of 15° — 30°. This value is, of course, 
much smaller than the maximum value of the angle 9i^2- 
This explains why MR(H) is larger for the interleaved 
configuration than for the separated configuration. 

The calculated results [Q for the separated configura- 
tion are given by the curves in Fig. 3. For each value of 
N, the three parameters ci^i{9max,i,i)'^, C2,2(^maa;,2,2)'', 
and Ci,2 were determined by fitting to the MR(H) data. 
The agreement between the calculated curves and the 
data is evident from the figure. 

To confirm that MR(H) for the separated configura- 
tion contains the contributions of [Co(10A)/Cu(200A)]n 
and of [Co(60A)/Cu(200A)]n, we also measured MR(H) 
for a multilayer containing only [Co(10A)/Cu(200A)]n 
and for another multilayer contain- 

ing only [Co(60A)/Cu(200A)]n. For each of these two 
multilayers, MR(H) consists of a single peak, located at 
the same magnetic field as one of the two peaks in the 
separated configuration. Thus, the two peaks observed 
for the separated configuration do indeed correspond to 
the two individual peaks. 

In conclusion, we have shown that the principal fea- 
tures of the MR(H) data can be explained quantitatively, 
for both the interleaved and the separated configurations, 
by invoking non-local electron scattering. 
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